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Abstract: In this paper we propose a new impact-parameter dependent CGC/saturation model. We 
introduce two new features in the model that make it consistent with what we know theoretically about 
the deep inelastic scattering. They are: the use of the exact form of the solution to the non-linear (BK) 
equation, whereas in all previous attempts only the form of r 2 Q 2 dependence, has been taken into account; 
and the large impact parameter dependence, through the 6-dependence of the saturation momentum which 
reproduce the correct behaviour of the amplitude at large impact parameters 6 {A oc exp (—//&)) as well as 
at large momentum transferred Qp (A decreases as a power of Qp as it follows from perturbative QCD). 
These improvement compared to all previous attempts to build such models, allows us to claim, that the 
experimental data are in accord with the prediction of CGC/saturation approach while previously, based 
on similar models, we could only conclude that the DIS data, perhaps, can be described by introducing the 
shadowing corrections at small photon virtualities. 
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1 Introduction. 

In this paper we continue our attempts to find an approach, based on Color Glass Condensate/saturation 
effective theory for high energy QCD (see Ref.[l] for the review), which will satisfy the following require¬ 
ments: it is simple but stems from the solution of the CGC/saturation equations. In other words, we 
would like to reduce the ad hoc or/and phenomenological ingredients to the minimum. We believe that the 
comparison of our model with a variety of experimental data will allow us to demonstrate how the collec¬ 
tive phenomena, that are incorporated in the CGC/saturation approach, are important for understanding 
scattering processes at high energies. 
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For comparison with experimental data, we choose the deep inelastic processes for which we have the 
CGC/saturation equations of Ref. [2] in a simple form of non-linear Balitsky-Kovchegov equation [3]. In 
addition to the simplicity of the theoretical interpretation of the deep inelastic structure function F 2 , these 
processes have been measured to high accuracy [4, 5], fitting to this observable allows us to determine all 
phenomenological parameters that are introduced in our model. However, the theoretical interpretation of 
this observable should also include the contribution of the colorless dipoles of large sizes, which require a 
non-perturbative description in the framework of QCD. In practice this means that we are not able to apply 
our approach to low Q 2 , where Q is the photon virtuality. It is advantageous to compare Fl with the data 
as the main contribution to this structure function, stems from dipoles of size 1/Q. However, only we have 
scant data [ 6 , 7] for Fl- 

After determining the values of all phenomenological parameters of our approach from comparing to 
the data on F 2 , we calculate and compare with the experimental data on open charm, and J/’F meson 
diffractive production data [8-11]. These data can be described within our theory with the same accuracy 
as F‘ 2 , as the wave function of J/T meson is known. However, other sources of uncertainty such as the 
influence of the real part of the amplitude and/or the corrections that stem from the skewedness effect, due 
to the fact that the gluons attached to the qq can carry different light-cone fractions x and x' of proton, 
have to be taken int account. They lead to uncertainties which are not under full theoretical control. The 
same remarks can be made regarding the virtual Compton reaction ( 7 * + p —>• 7 + p) [12] for which the 
wave function is known to even better accuracy than for J/'F production. 

From the theoretical point of view, the description of the diffractive production of 4> and P resonances 
[13] is less reliable, since we do not have a theoretical approach for the wave functions of these mesons, this 
in addition to all other uncertainties. 


The main goal of comparing with experimental data, is to show that the CGC/saturation approach is 
able to describe all these processes. This paper is neither the first nor the last effort to build such a model 
based on the equations of CGC/saturation approach. Numerous attempts have been undertaken over the 
past two decades (see Refs. [14-33]) to build such models. Therefore, we clarify, in the introduction, the 
aspects of our model which are different. 

First of all, it is necessary to construct a model, since the CGC/saturation equations cannot reproduce 
the correct behavior of the scattering amplitude at large impact parameters (see Ref. [34, 35]). Such failure 
leads, at least, to two conclusions: first, we cannot trust the solution of the CGC/saturation equations, and 
second, in attempting to describe experimental data we are doomed to build models to introduce the main 
features of the CGC/saturation approach. 


We introduce the non-perturbative impact parameter behavior in the saturation moment, accordingly 
to the spirit of geometric scaling behavior of the scattering amplitude[36] and according to the semi-classical 
solution to the CGC/saturation equations [16]. Similar assumptions for the non-perturbative 6 -behavior of 
the scattering amplitude, is typical most models on the market (see Refs. [21-25, 28, 33]) and we refer only 

to the behavior at large 6 , of the saturation scale Q s (x, b). We introduce Q s (x, 6 ) exp (—pb) while 

in all other models Q s (x,b) — exp (—/r 2 6 2 ). It should be stressed that the exponential decrease at 
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large b, is in accord the Froissart theorem [37]. We would like to emphasize that the b dependence in our 
model reproduces the perturbative QCD prediction for the large values of the momentum transferred (Qt) 
leading to power-like behaviour of the scattering amplitude at large Qt, which cannot be reproduced for 
the Gaussian b distribution in other models. 

The main difference stems from the way we find the behavior of the dipole-nucleon scattering amplitude 
in the saturation region. In general, we follow the procedure suggested in Ref. [18]: we assume the geometric 
scaling behavior of the amplitude in the saturation region, and match the asymptotic behavior [38] deep 
inside the saturation region, with the behavior of this amplitude in the vicinity of the saturation scale 
[1, 39, 41]. We also take into account, not only the correct behavior of the scattering amplitude inside 
the saturation region (see Ref.[38]): viz. N oc exp (—Az 2 ) with z = In ( r 2 Q 2 ) where r is the size of the 
dipole and Q s is the saturation momentum, but the correct value of the coefficient A [38] and we specify 
the procedure of matching at z = 0. The procedure will be described below, and it is based on Ref. [42], to 
which we refer the reader for all details of this approach. 

The natural question arises whether we introduce only slight, cosmetic changes to the models of 
Refs. [21-25, 28, 33] or the changes are principal in nature. We introduced in the model all theoretical 
information about the deep inelastic scattering, as, we believe, should be done in all phenomenological 
models, since the available experimental data are not sufficient to allow us to differentiate between theo¬ 
retically correct and theoretically insufficient (or even wrong) approaches. It is enough to give an example 
of the deep inelastic scattering which can be described without any shadowing corrections (e.g. Ref. [40]). 
Therefore, we need only demonstrate that the theoretical changes that we introduce, are large enough to 
influence the description of the experimental data. This is the case since the correct 5-behaviour leads to 
the ln 2 (l/x) increase of the deep inelastic structure function F 2 ,while the correct behaviour at large z 
determines the coefficient in front of ln(l/x) in F 2 . Our numerical estimates show that both factors are 
essential in the description of the HERA data, and become even more significant at higher energies. 

As will be seen from the body of the paper, our model provides a good description of the experimental 
data. The fact that this description does not appear better than in the other insufficient models on the 
market, allows us to claim that the CGC/saturation approach does not contradict the available experimental 
data. On the other hand we learned only that shadowing corrections are, perhaps, needed to describe the 
experimental data from comparison with the experiment all models of Refs. [14-33]. 

2 Theoretical input 
2.1 General formula 

The general formula for deep inelastic processes takes the form (see Fig. 1 and Ref. [1] for the review and 
references therein) 

N (Q, Y ; b) = J jT 1 dz %. (Q, r, z) N (r, Y ; b) ** v (r, z) (2.1) 
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where Y = In (1 / xBj ) and XBj is the Bjorken x. z is the fraction of energy carried by quark. Q is the 
photon virtuality, b is the impact parameter of the scattering amplitude. 

Eq. (2.1) splits the calculation of the scattering amplitude into two stages: calculation of the wave 
functions, and estimates of the dipole scattering amplitude. 



JAU, cj>, p) 


Figure 1. The graphic representation of Eq. (2.1) for the scattering amplitude, 
of the interacting dipole, z is the fraction of energy that is carried by one quark, 
scattering amplitude 


Y = ln(l /xsj) and r is the size 
b is the impact parameter of the 


2.2 Dipole-nucleon scattering amplitude 

We follow the general plan of Ref.[18], searching for a solution to the CGC/saturation equation, by matching 
the solution in the vicinity of the saturation scale, with the solution deep inside the saturation domain. In 
both these regions we have analytical solutions. 


2.2.1 Vicinity of the saturation scale 

Dipoles of sizes r are such, that r 2 Q 2 (Y, b) —>• 1 the scattering [1, 39, 41] 

N(r,Y;b ) = N 0 (r 2 Q] (h)) 1_7or 
where 7 cr is the solution of the equation 

X {lev) = dx (7 cr) 

1 7 cr d^cr 

In Eq. (2.3) x ( 7 ) is the BFKL kernel that takes the form 

xin) = 2^(1) - V’(t) - V’C 1 - 7) 


( 2 . 2 ) 


(2.3) 


(2.4) 


where ^ (z) is the digamma function. 

The scattering amplitude of Eq. (2.3) shows the geometric scaling behaviour [36]. In Ref. [39] it is 
shown that the first corrections due to a violation of the geometric scaling behavior can be taken into 
account by replacing 1 — 'y cr in Eq. (2.2) by the following expression 

1 - 7cr -> 1 - 7cr - 2k \y ln ( r2 ( 6 )) ( 2 - 5 ) 

where A = a s (x (7cr) / (1 - 7cr)) and k = x" her) /x! her)- 
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2.2.2 Saturation region 


The new ingredient in our model is the way we found the solution to the Balitsky-Kovchegov (BK) equation 
[3] in the saturation region where r 2 Q 2 ( Y,b ) > 1. In this region the scattering amplitude demonstrates 
geometric scaling behavior, being a function of one variable 

* = In (r 2 Q 2 (T,6)) = In (r 2 Q 2 (Y = Y 0 ,b)) + A (Y - Y 0 ) = £ + A (Y — Y 0 ) (2.6) 

since the saturation scale behaves as 

Ql (Y, b) = Q 2 (Y = Y 0 ,b) e x {Y ~ Yo) = Q 2 (Y = Y 0 , b) (^)* (2.7) 


where Y = ln(l/x), x is the Bjorken energy variable which is equal to x = Q 2 /s for deep inelastic scattering. 
Q is the photon virtuality and s is the energy squared for the process. 

Deep in the saturation region (z l,r 2 Q 2 (Y,b) 2> 1) the solution to the non-linear BK equation 
is known (see Ref. [38]). In Refs. [42] the next order correction to asymptotic behavior at large z has been 
calculated. In this paper we use the following expression taken from Ref. [42] 


iV 2>1 (z) 

Z 


l-27le- 2 - V2XA 2 ~^e~ 2Z 

(z - \A\J A7t/2 - 2^(1)) 

2A 


+ 0(e ~ 3Z ) 


( 2 . 8 ) 


where tp(x) is the digamma function (see Ref. [43] formula 8.360 - 8.367). 

The second term in Eq. (2.8) is the solution given in Ref. [38] in which the theoretically unknown 
constant A is introduced both as the coefficient in front and as correction to the argument. The third term 
is the next order correction at large z. 

For z > 0 but z <C 1, the solution to BK equation has been found in Ref. [42] and it takes the form 


iV' 


0 <z«l 


(z) = -2 


A (1 - 7cr ) 2 (IVo - A (1 - 7cr )) N 0 e-a-7^)* 


{N 0 - (N 0 - A(l- 7cr ))e-( 1 -7c r)zy 
where A^o is the value of the scattering amplitude at z = 0. 

We match these two solution at z = z m where 

dAr o<z«i ( z = Zm) dN z^i (z = Zm ) 


(2.9) 


N 


0<z<l 


(z = Zm) = AI 2>i (z = Z m ) \ 


dZri 


dz„ 


( 2 . 10 ) 


It has been demonstrated in Ref. [42] that we find z m ~ 1. 
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2.2.3 Estimates for the parameters of the scattering amplitude 


The entries in Eq. (2.2), Eq. (2.8) and Eq. (2.9) can be divided in two groups: parameters that can be 
evaluated in the leading log(l/x) approximation of perturbative QCD, and phenomenological parameters 
that can be found only by comparing with the experimental data. The value of 7 cr is determined by Eq. (2.3) 
and it is equal to 0.37. The energy dependence of the saturation scale is given by 

A = dsx(7cr)/(l-7cr)- (2.11) 

In Eq. (2.5) k = x" (7 cr) / x' (7 cr) in the leading order of perturbative QCD(LO). In our procedure of fitting 
we consider the value of A as a parameter of the fit since we do not know the value of as and the next-to- 
leading order corrections (NLO) in perturbative QCD change considerably the value of A [41, 44] making it 
in 3-4 times less than the LO estimates. The value of A is determined by the linear BFKL equation which 
is known in the NLO, and therefore, we consider the estimates of A as reliable, which have to be taken into 
account in spite of uncertainties in the value. A ~ 0.2 which we obtain from the fit is close to the NLO 
estimates. 

In principle the value of Nq can be calculated using the linear evolution equation with the initial 
conditions. However, it depends on the phenomenological parameters of this initial condition. We choose 
to use Nq as the parameter of the fit. 

The impact parameter dependence of the saturation scale as we have discussed, is the pure phenomeno¬ 
logical input, which can only be estimated in non-perturbative QCD. For Q s (Y = Yq, b) we use the following 
expression 

Q 2 (Y = Y 0 ,b ) = m 2 S (b) = m 2 (mbK 1 (mb)) 1/{1 -' ycr) (2.12) 

m must be determined from the fit. We differ from other models in that Eq. (2.12) leads to 
Q 2 (Y = Yq, b) exp (— mb/(l — 7 cr )) providing the correct large b behavior of the scattering ampli¬ 

tude. In other models (see Refs. [21-25, 28, 33]) Q 2 (Y = Yo,b) oc exp (—b 2 /B). 

2.3 Wave functions 

The wave function in the master equation (see Eq. (2.1)) is the main source of theoretical uncertainties: 
even in the case of deep inelastic processes, we can trust the wave function of perturbative QCD only, 
at rather large values of Q 2 > Qq with Q q ss 0.7 GeV 2 (see Ref. [45]). The expression for ('L*'L) 7 * = 
T 7 * (Q,r, z) 'L 7 * ( Q,r,z ) is well known (see Ref.[l] and references therein) 

(T*^)J = {[z 2 + (1 - z) 2 ] e 2 Kf(er) + m 2 /to(er)} , (2.13) 

(***)? = —a em Y j e)Q 2 z 2 (l-z) 2 K 2 (er), (2.14) 

f 

where T(L) denotes the polarization of the photon and / is the flavours of the quarks, e 2 = + Q 2 z( 1 — 

*)■ 
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In addition to the total DIS cross section, the wave function is known theoretically for deeply virtual 
Compton scattering (DVCS). For this process ^*^/) DVCS = (r,z) T 7 * ( Q,r,z ) is equal to 

o n 

('F*'F)y' c ' s = { [z 2 + (1 - z) 2 } eiVi (er)m/iFi (m/r) + m 2 K 0 {er)K 0 (mfr)} . (2.15) 


Since in DVCS the real photon is produced, so its polarization can be only transverse, as seen from Eq. (2.15). 

For vector meson diffractive production, the wave functions are known only for the mesons that are 
constituted of heavy quarks. The most popular example is the J/\F diffractive production, but we have to 
bear in mind, that the mass of the charm quark is not very large, and corrections can be essential. For all 
other mesons the confinement corrections are large and the form of the wave functions motivated by the 
heavy quark mesons, can be considered only as a pure phenomenological conjecture. We use the following 
form of the overlaps integrals (' Vy^)T,L = 1F 7 * (Q, r ,z) (r, z) taking them from Ref.[22] 


('FO'F)t = e/e ——- {m 2 fK 0 (er)^T(r, z) - [z 2 + (1 - z) 2 ] eiti(er)<9 r </V(r, z)} , 
Tizyi — z) J 


(Vv'P)l = e/e — 2Qz(l - z) K 0 (er) 


7r 


mn 2 f — V 2 

M v (j) L (r , z) + ■ - t </> l(d z) 


Myz(l — z) 


m 2 f n 2 2z(l - z)r 2 rnjTt 


2 - 772 ' 


<t>T,L{r , z) = Mt.lz(\ - z) exp - 


8*(1 - z) 


n 2 


+ 


(2.16) 

(2.17) 

(2.18) 


where the effective charge e/ = 2/3, 1/3, or l/\/2, for J/t/, q i, or p mesons respectively, rrif = 140MeV is 
used in all above equations and the values of Mt,l and 1Z are taken from Table 2 in Ref. [19]. 


2.4 Other theoretical uncertainties 


2.4.1 Real part of the scattering amplitude 


The master equation (see Eq. (2.1)) is derived assuming that the scattering amplitude is purely imaginary. 
To account for the real part of the amplitude, we need to multiply the differential cross section calculated 
by using Eq. (2.1) by factor 1 + p 2 where p is the ratio of real to imaginary parts of the scattering amplitude 
A. In principle, we can calculate the real part of the amplitude if we know the imaginary part using the 
dispersion relation. For large values of energy this calculation simplifies to 


7 T din (IV) 
2 din(1/x) 


(2.19) 


for a more general formula for the amplitudes that show Regge-type behavior at high energies (see Ref. [46] 
and references therein) 

8^[Npr^) 

<91n(l/x) 


p = tan(7r/3/2), with /3 = 


( 2 . 20 ) 












2.4.2 Off diagonal contributions 


For vector meson production or DVCS, we need to modify the main formula of Eq. (2.1), to account for 
the fact that the virtual photon, and the produced vector mesons carry different fractions x and x' of the 
proton’s (light-cone) momentum. In the leading ln(l/x) limit, this imbalance (skewed effect) does not affect 
the scattering amplitude, but it has to be taken into account in the limit that x' <C x <C 1. Two include 
the skewed effect we multiply the amplitude by the factor Rq given by [47] 


where /3 is given by Eq. (2.20). 


RgW) 


2 2 / 3+ 3 r (/3 -f 5/2) 

r(£ + 4) 


( 2 . 21 ) 


3 Fitting F 2 and values of the parameters 

The most accurate experimental data available are on the deep inelastic structure function F 2 [4]. Our 
strategy is to determine all phenomenological parameters by fitting to these data. Other data we will 
use to illustrate that our model can be used as a simple and convenient tool in describing a variety of 
processes. As has been mentioned, we can trust our model in the restricted kinematic region, which we 
choose in the following way: 0.85 GeV 2 < Q 2 < 60 GeV 2 and x < 0.01. The lower limit of Q 2 stems from 
non-perturbative correction to the wave function of the virtual photon, while the upper limit is originated 
from the restriction x < 0.01. 


A 

N 0 

m {GeV) 

Ql {GeV 2 ) 

m u { MeV) 

mrf(MeV) 

m s (MeV) 

m c (GeV) 

X 2 /d.o.f. 

0.197 

0.34 

0.75 

0.145 

2.3 

4.8 

95 

1.4 

178/155 -1.15 

0.184 

0.46 

0.75 

0.118 

140 

140 

140 

1.4 

176/154 = 1.14 


Table 1 . Fitted parameters of the model. Qq = m 2 Xq. 


We have the following fitting parameters: A for ^-dependence of the saturation momentum, Nq for the 
value of the scattering amplitude at r 2 Q 2 = 1, m for the impact dependence of the saturation scale and 
xq or Q q = 171 2 Xq for the value of the saturation scale. Masses of the quark we do not regard as fitting 
parameters and consider two sets of these masses. In the first set we take the current masses (see the first 
row of Table 1) and we consider this as the most reliable fit based on the consistent theoretical approach.We 
also make a fit putting all masses of light quarks (second row of Table 1) to be equal to 140 MeV. We view 
this mass as a typical infra-red cutoff , which was used in the wave functions of the produced mesons ( rrif — 
140 MeV in Eq. (2.16) and Eq. (2.17)). 

Table 1 gives the values of the fitting parameters and Fig. 2 demonstrates the quality of the fit. 

One can see from the Table 1 that both sets of parameters give good descriptions of the experimental 
data and the values of the parameters are very close for both fits. The value of A leads to very small value 
of the QCD coupling as in Eq. (2.11), but, as we have discussed, the estimates in the NLO of perturbative 



















X 

Figure 2. Our fit of F 2 with the values of parameters given by Table 1. The first set of parameters is shown in solid 
red curves while the second in blue dotted lines. The data is taken from Ref.[4], 


QCD lead the value of A[41, 44] which is much smaller than in Eq. (2.11) and which is close to the value 
that we obtain from the fit. The value of Nq turns out to be small enough, to use our matching procedure. 
Qq = m 2 Xq looks reasonable and generates the values of the saturation momentum shown in Fig. 3. 

The value of m is smaller than the typical mass in the electro-magnetic form factor of the proton, 
but we do not expect that it will be the same. The mass that we extract from our fit is close to the 
mass of p-meson. However, it should be noted that the decrease of Q 2 at large b is proportional to 
exp (—j^—b'j = exp ( —1.2 (GeH -1 ) b ). On the other hand the behavior of amplitude with b differs 
from the saturation scale. In Fig. 4 one can see that both the saturation, and the violation of the geometric 
scaling behavior influence the resulting b-dependence of scattering amplitude. The saturation flattens 
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Qfa, b ) 



- b=0 

- b=0.2 fm 

- b= 0.4 fm 

b=0.7 fm 
- b= 1 fm 


Figure 3. The value of the saturation momentum Q 2 (x, b) versus x at fixed b for the parameters given by Table 1. 


the 6-dependence at small values of 6, while the large b behaviour shows a more rapid decrease than the 
6-dependence of the saturation scale (see Fig. 4). 



Figure 4. The 6-dependence of the scattering amplitude for the parameters given by Table 1. S (b) is given by 
Eq. (2.12). 


It should be stressed that in the framework of our parametrization of the 6-dependence of the saturation 
momentum, the scattering amplitude decreases as exp (—mb) while in all other models on the market it has 
a Gaussian behavior: exp (—m 2 6 2 ). 

Fig. 5 we present the comparison between our fit of F 2 with two sets of parameters at low values of Q. 
The set with large masses of quarks leads to a much better description illustrating the the non-perturbative 
corrections to the wave function of the virtual photon are essential at Q 2 < 0.85 GeV 2 . 
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Figure 5. The ^-dependence of F£ c at small values of Q 2 < 0.85 GeV 2 for the parameters given by Table 1. The 
red (upper) line corresponds to set 1 (upper row of Table 1) while the blue one (low) is the description with set 2. 
The data are taken from Ref. [6, 7]. 

4 Comparison with the experimental data 

In this section we compare the predictions of our approach using the set of parameters shown in Table 1, 
with the experimental data. 

F^ 0 : The contribution of the cc pair to the deep inelastic structure function can be calculated with 
the same theoretical accuracy as the inclusive F-i. In Fig. 6 we compare the HERA data on F£° [5] with 
the theoretical predictions. One can see that the agreement is reasonable. 


Fl ■ Fl can be calculated within the same accuracy as F% 5 , and the comparison with the scant data 
available [6, 7] is plotted in Fig. 7. Two sets produce the same quality of the descriptions since the values 
of Q are rather large. 

Main formulae for cross sections of exclusive diffractive production: 

First, we need first to calculate the scattering amplitude at fixed t which takes the following form [22]: 


A 


7*p—>Vp 
T,L 


(x, Q, A) = i j d 2 r J^J d 2 b )t,l e^-Md-A N (r> y. h) (4.1) 


where |A | 2 = —t. 
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10 4 10 2 


X 

Figure 6. The ^-dependence of Fg c at fixed values of Q 2 : 0.85 < Q 2 < 60GeE 2 for the parameters given in Table 
1. The data are taken from Ref. [5]. 


The elastic diffractive cross section is then given by 


d 7 *p-+Vp 
aa T,L _ . 7 * p ^y p 2 

d t 167T T,L 


1 

16-7T 


f d2r f 0 ^ J d2b N (r, Y ; b) 


. (4.2) 


If we restrict ourselves to measuring two integrated observables: 

7 *p->Vp \ 

-1t=o ), Eq. (4.1) and Eq. (4.2) can be simplified. 


] 7 *P^Vp 

= f dt T, dt - and Bn 


In ' drTT ’ L 


dt 
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Q 2 =3|.5GeV 2 

-A~ 

Q 2 =12G 

■ ij 

eV 2 

i 

'■- " 

Q 2 =20GeV 2 

""f f— : 

Q 2 =45GeV 2 

: 

Q 2 =60GeV 2 
; -! r ; 

10 3 

F l 

• HI 

• ZEUS 


10 3 10 3 
X 


Figure 7. The x-dependence of Fp at fixed values of Q 2 : 0.85 < Q 2 < 60 GeV 2 for the parameters given in Table 
1. The red (blue) lines correspond to set 1 and set 2 fits. The data are taken from Ref. [7]. 


The formulae for them take the following form 


a 


ry*p—>Vp _ 


1 

167T 


d 2 b 


(4.3) 


{ J d2r Jo^^ )T ' L N ^ Y ^}{ j dV iV(r',y;b-(l-/y)}^ 


a 


*p^ Vp B d = 


167r 


d 2 6d 2 r d 2 r' 


r x dz r l dz' 

lo 4tt 7 0 4tt 


(4.4) 


(46 2 + (1 - zfr 2 + (1 - z') 2 r' 2 ) {(^) t ,l N (r, T; b) } { (4^T)t,l N (/, Y;b — (1 — z)r') 
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Deeply virtual compton scattering(DVCS) 

Deeply inelastic Compton scattering: 7 * + p —>• 7 + p, can be calculated to the same accuracy as all 
reactions that has been discussed, but it suffers from the errors due to the procedure for calculating both 
the real part of the amplitude and its skewness. Fig. 8 shows the comparison with the data, taken from 
Refs. [12]. One can see that the description of the cross sections are satisfactory, while the energy dependence 
of the slope turns out to be more complicated than in the experimental data. 


* 



Y*P->YP Y*P->YP 



Fig. 8 -a 


Fig. 8 -b 


Fig. 8 -c 


Figure 8. 

scattering, 
to the set 1 


The energy (Fig. 8 -a) and Q 2 (Fig. 8 -b) dependence of the cross section for deep inelastic compton 
Fig. 8 -c shows the energy dependence of the slope. The data are taken from [12] . The lines correspond 
of the parameters. The red color describe the HI data while the blue one stand for ZEUS data. 


Total cross sections for diffractive production of vector mesons a 1 p ^ Vp (see Eq. (4.3)) 

I 11 Fig. 9 and Fig. 10 we compare our prediction with the experimental data [8-11] for productions of 
J/\F, (j) and p-mesons. 

We introduce the factor (1 + P 2 )Rq to account for the real part of the amplitude and for the skewness 
as has been discussed in section 2.4. One can see that agreement is good, and the heavier the produced 
meson, the better the agreement, as expected. 

The slopes Bp: F>d are calculated using Eq. (4.4) and Fig. 11 shows the comparison with the 
experiment data taken from Refs. [8-11]. In Fig. 12 we plot the dependence of the slope on the energy 
W. We are able to describe the values and main regularities of the slope’s behavior on Q 2 and W. Fig. 12 
shows the the shrinkage of the diffraction peak is very mild. 
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y p — > p p 


Y* p -»J/v|/ p 



y* p -> <|) p 




Fig. 9-a 


Fig. 9-b 


Fig. 9-c 


Figure 9. The VF-dependence of a 1 p^ v p for vector meson production at fixed values of Q 2 . The data are taken 
from [8-11]. 


y* p —> J/xj/ p YP^^P y* P P P 





Figure 10. The Q 2 + M^-dependence of er 7 p^ v p for vector meson production at fixed values of W. The data are 
taken from Refs. [8-11]. 


5 Conclusions 

In this paper we add to the many CGC/saturation models [14-33] two crucial ingredients: the correct 
solution to the non-linear (BK) equation [3] in the saturation region, and impact parameter distribution 
that leads to exponential decrease of the saturation momentum at large impact parameters and to power- 
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y* p -» JA|/ p y* p -»<|) p y* p ^ p p 





Fig. 11-a Fig. 11-b Fig. 11-c 

Figure 11. The Q 2 -dependence of Bp for vector meson production at fixed values of W. The data are taken from 
Refs. [8-11]. 


Y* P -> P P 



Figure 12. The TF-dependence of Bp for p production at fixed values of Q 2 . The data are taken from Refs. [8-11]. 


like decrease at large transferred momentum that follows from perturbative QCD. We obtain the solution to 
the BK equation based on the ideas proposed in Ref. [18]: to match the analytical solution for the scattering 
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amplitude deep inside of the saturation domain [38] with the solution in the vicinity of the saturation 
momentum [39, 41]. The novel features of our approach is that we use the exact form of the solution 
to the BK equation (see section 2.2 and Refs. [38, 42]) but not the form of r 2 Q 2 dependence as have bee 
assumed in all previous attempts. We introduce the large impact parameter dependence through the in¬ 
dependence of the saturation momentum. The difference between our and the other attempts consists in 
the assumption that Q 2 oc exp (— mb ) at large b, while in all other models the Gaussian behavior at large 
b Q 2 oc exp (— m 2 b 2 ) is assumed. In the framework of our model, the exponential b dependence of Q 2 
translates in the exponential decrease of the scattering amplitude at large b, satisfying the Froissart theorem 

[ 371 - 

Using four fitting parameters we obtain good overall agreement with the experimental data. Comparison 
with the experimental data we found two regularities. The first one is , the better we know theoretically 
the wave functions of the produced mesons, the better the agreement with experimental data. The second 
is the fact, that the energy dependence of the saturation momentum, is much milder that it is predicted by 
the BK equation in leading log(l/x) approximation of perturbative QCD. The last observation, we hope, 
will stimulate the search for the non-linear equation, in the next-to-leading order (NLO) of perturbative 
QCD. The first attempts to estimate the energy behavior of the saturation scale in NLO show that the 
value of A significantly decreases (see Refs. [41, 44, 49, 50]). 

We believe that the model presented here, will be a useful tool to estimate the CGC/saturation effects 
in the variety of exclusive reactions. 

The model includes everything that we know from theory about deep inelastic processes and, as such, 
can be used for the comparison at higher energies including the LHC energies. From a good description 
of the experimental data, which is not better than in insufficient models of Refs. [14-33] we can conclude 
that the CGC/saturation approach does not contradict available experimental data. The weak statement 
is nevertheless stronger than the one that we can make from the models of Refs. [14-33]: shadowing 
corrections are, perhaps, needed to describe the experimental data. Formulating the result of the paper in 
one sentence, we firmly believe that only our model is reliable for extrapolation to higher energy including 
the LHC one. 

We need to discuss why we use a model while, at first sight, we have a good description of the experimen¬ 
tal data at least for the deep inelastic scattering, based on the leading order(LO) or even beyond of the lead¬ 
ing order (NLO) perturbative QCD evolution equations (see Refs. [50-52]). Unfortunately, these equations 
as we have discussed, cannot provide reliable predictions for physical observables including F 2 , while the 
problem of large impact parameter behaviour of the scattering amplitude would not be solved theoretically 
and the CGC/saturation equations would have to be be modified to include this behaviour. The solutions 
of the CGC/saturation equations that have been discussed, are valid only “ under the approximation that 
the dipole scattering amplitude is independent on the impact parameter" (citation from Ref. [52]). Such an 
approximation generates the dipole scattering amplitude N oc 6^ (Qt) and F 2 oc ImiV (Qt = 0) —>• 00 , the 
prediction which cannot be compared with the experimental data. Therefore, we doomed to build models 
for comparing with the experimental data. Using the results of these models we can evaluate how essential 
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is the 6-dependence for the description of the experimental data, and see that it is essential even in the 
HERA kinematic region, and will be even more important at higher energy (i.e. at the LHC). 

The second question is why we improve the model based on the leading order perturbative QCD while 
we know that the next-to-leading QCD corrections (NLO) are essential. The answer is: first we need to 
build a reliable model in the leading order QCD and only after this we can build the model that takes into 
account the NLO corrections. We believe that the NLO corrections will change all parameters of the model 
especially the energy behaviour of the saturation scale leading to small values of A as it was demonstrated 
in Refs. [41, 49, 50]. Since the value of A can be found from the solution of the linear equation we trust 
these estimates and, as we have discussed, we see in this NLO result the support for the phenomenological 
value of A. However, we believe that we do not have enough solid theory information to discuss in the 
NLO the change in the behaviour of the scattering amplitude inside the saturation region building a model 
that includes the impact parameter behaviour. The general form of the non-linear equation in the NLO [53] 
shows that we can use the approach of Ref. [38] for finding an equation at * > 1 which takes the form 

— = -as J d 2 z^K LO (x,y\z) +a s K NLO (x,y\z) ^N(Y;x,y) (5.1) 

K nl o t urns out to be such that K LO (x,y\z) + as K NLO (x,y\z) —> 0 for asz ~ 1[44, 49, 50]. It means 
that we cannot specify large z asymptotic behaviour using the NLO kernel. We have to re-sum all orders 
of ( asz) n but we do not know a theory proof of the form of the non-linear equation for such a re-summed 
kernel. All attempts to write the non-linear evolution equation for the re-summed kernel (e.g. Refs. [49, 50]) 
use the simplified form of the kernel for the BFKL equation taking into account only its leading twist 
contribution. The correct coefficient from the full BFKL kernel in front of z 2 term in the solution at large 
z 1, is much more important for phenomenology than the NLO corrections. Recalling that the running 
QCD coupling that should be included in such a re-summation, leads to the violation of the geometric 
scaling behaviour [54], and that the improved double log approximation of Ref. [50] cannot be used in the 
saturation region[38, 44, 49], we see that a lot of theory work must be done before we will be able to build 
a model that will include the impact parameter behaviour and the NLO corrections to the BK equation. 
We plan building such model as a natural goal for future work. 
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